Rabi Oscillations in Landau Quantized Graphene 
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The canonical model of quantum optics, the Jaynes-Cummings Hamiltonian describes a two level 
atom in a cavity interacting with electromagnetic field. Graphene, a condensed matter system, 
possesses low energy excitations obeying to the Dirac equation, and mimics the physics of quantum 
electrodynamics. These two seemingly unrelated fields turn out to be closely related to each other. 
We demonstrate that Rabi oscillations, corresponding to the excitations of the atom in the former 
case are observable in the optical response of the latter in quantizing magnetic field, providing 
us with a transparent picture about the structure of optical transitions in graphene. While the 
longitudinal conductivity reveals chaotic Rabi oscillations, the Hall component measures coherent 
ones. This opens up the exciting possibility of investigating a microscopic model of a few quantum 
objects in a macroscopic experiment of a bulk material with tunable parameters. 

PACS numbers: 

Graphene, a single sheet of carbon atoms in a honeycomb lattice has attracted enormous interest recentljii. Its 
quasiparticle states are pseudospinors, where the spinor structure is a consequence of the two-fold sublattice structure 
of the honeycomb lattice (Fig. [1]). They obey a two-dimensional Dirac equation, whose speed of light is replaced 
by the Fermi velocity (being l/300th the speed of light). This implies a number of striking properties, including the 
unconventional quantum Hall effecli^, Klein tunneling^ and Zitterbewegung^i^ due to particle-hole excitations. In 
an applied magnetic field, perpendicular to the carbon sheet, the formation of Landau levels En with an unusual 
dependence En ^ y/ri on the Landau level index n was predicted, and also observed experimentallj*2i^. 




FIG. 1: (Color online) A small fragment of the honeycomb lattice is shown, the filled red and empty black circles denote the 
two sublattices. 

Coupling of a (pseudo)spinor to an external quantum field is a common feature in quantum optics, which deals with 
the interaction of electromagnetic field and matter. In this context, one of the simplest, fully quantum mechanical 
model is the Jaynes-Cummings (JC) Hamiltonian^ 

Hjc = V{a+<T- + (T+a) + Acr, , (1) 

which describes a single two-state atom, represented by the Pauli matrices, that is interacting with a (single- mode 
quantized) electromagnetic fieldi^. a"*" (a) are the photon creation (annihilation) operators, and V is the coupling 
strength between the atom and the electromagnetic field. The interaction leads to a periodic exchange of energy 
between the electromagnetic field and the two-level system, known as Rabi oscillations. This effect can be interpreted 
as a periodic change between absorption and stimulated emission of photons. It is simply a consequence of the fact 
that the atom is not in an eigenstate after the absorption of n photons. A very similar situation occurs for the Landau 
levels of graphene. These levels play the role of the atomic levels in quantum optics due to photon absorption and 
the sublattice structure of graphene provides the two-fold degeneracy of the two levels of the atom. In terms of the 
Hamiltonian ([T|), the photon operators a+ and a are operators that are acting on Landau levelsii: a+ = 7r+/ \/2eB 
and a = Tr~ /\/2eB, where tt^ = ttx i TTy with tt = p -|- eA are the (Peicrls substituted) momentum operator in the 
vector potential A. Formally, they act as creation and annihilation operators of a harmonic oscillator and obey a 
bosonic commutation relation [7r~,7r+] = 2eB. Moreover, the coupling constant V depends on the magnetic field B 
and the Fermi velocity in graphene vp = 10^ m/s through V = vp\/2eB. Finally, A represents a possible excitonic 
gap>i^ or a substrate induced bandgajJ^ in epitaxial graphene. 
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The pseudospinor state in the nth Landau level is not an eigenstate of the quasiparticle Hamiltonian of graphene. 
If we prepare an initial quasiparticle state of graphene in a certain Landau level by switching on an external magnetic 
field, the quasiparticles will go from the initial non-equilibrium state to other Landau levels. Depending on the 
damping due to (impurity) scattering, the quasiparticles oscillate between different Landau levels until they reach 
their equilibrium state, which is a superposition of Landau levels. Recent experiments on high-mobility samples of 
graphene have indicated that damping due to impurity scattering can be negligible on experimentally relevant time 
scalesH*i5. 

In this Letter we shall discuss the non-equilibrium current dynamics of high- mobility graphene. Quantum states are 
prepared by a short electric pulse and the subsequent current dynamics is controlled by oscillations between Landau 
levels. 

The Hamiltonian Hjc can be diagonalized for Landau levels, and the resulting eigenvalues are obtained as Ena = 
ay^ A'^ + V^{n + 1), where n = 0, 1, 2. .. non-negative integer, a ~ ±. In the non-relativistic limit (A 3> V), 
the usual Landau level spectrum is obtained as a(A -I- LUciji + 1)) with the cyclotron frequency of massive (Dirac) 
fermions uJc = VpcB/A. In addition, there is a special eigenstate, stemming from the Landau level at the Dirac 
point with E* ~ —A, which formally corresponds to n = —1 and a = —1. Having determined the spectrum of 
the Hamiltonian, we turn to the investigation of the current correlations. Since cr^f^y-^ is the current density in x 
(jj) direction due to the equation of motion jx{y) = i[H,x{y)] ~ VFO'x{y)j we can calculate the dynamical correlation 
function Cxx{t) ~ {o'x{t)ux{0)) and Cxy = {<7x{t)(Jy{0)) (symmetric and antisymmetric dipolc-dipole correlator). In 
quantum optics, these describe the transitions between the two atomic states, and tell us about the spectrum of 
Rabi oscillationsi^. On the other hand, Cxx{t) and Cxy{t) play the role of the longitudinal and Hall current-current 
correlation functions in the Dirac case, and leads eventually to the optical conductivity of graphen o^^'^^ . Therefore, we 
expect the well-known Rabi oscillations of quantum optics characterizing the excitations of the atom to be observable 
in the response functions of Landau quantized Dirac fermions^. 

We start with the general correlator C{(p,t) = {ax{t){cos{(p)ax + s'm{ip)ay)) , evaluated using the (structureless) 
bosonic a operators. This defines both Cxx{t) — C{Q,t) and Cxy(t) = C(7r/2,t), and reads as 

C(V3, t) = ^ gcfiEna) CXp{i{Ena + En-sf)t + is(p)Pna-*n~sf + QcfiE*) ^ exp(i(£'o7 + E*)t - iip)Pt->0'r, 

(2) 

where f{E) = l/(exp((i? — l^)/T) + 1) is the Fermi function, gc = NfAceB/2TT is the degeneracy of the Landau levels 
and spins, Ac is the area of the unit cell, to be taken as unity in the Dirac approach, Nf = 2 stands for the spin 
degeneracy. From this, it follows immediately that a E^a Landau level with n > and given a possesses 4 possible 
optical transitions to the adjacent levels as En±i±a (on the same side and on the other side of the Dirac cone), the 
n = level 3 transitions to Ei±a and E* and the E* level two transitions to Eo±. The non-zero transition matrix 
elements are given by 

Pnc^n-s-, = \{^ + 1^) (l + ;^^) f"'^ " ^ 0, and P,^07 ^\{^- ' (3) 

which satisfy X]m7 Pna^m-y ~ 1, arid agree with the transition probabilities for Rabi oscillations of atoms induced by 
external electromagnetic field. These approach 1/4 in the classical limit (of bosons) n — > oo, in which case the field 
contains many bosons, whose quantum character can then be neglected^. Interestingly, the A = limit yields the 
classical matrix elements for any n ^ 0. However, the E* level never reaches the classical limit, and is responsible for 
the anomalous optical properties of graphene in magnetic field^i. The evaluation of the microwave conductivities of 
graphene follows readily from Eq. and the response function is Xxa{t) = ~Q{t)2e^VplmCxa{t) , which yields to 
the optical (a = x) and microwave Hall (a = y) conductivity as axaiijj) = Xa:a(w)/«a; upon Fourier transformation. 

Eq. ([2]), together with the relation to the JC Hamiltonian provides us with a particularly simple picture about the 
optical selection rules and transitions by relating them to the Rabi oscillations. Therefore, the optical conductivity by 
varying the frequency sweeps through all possible transitions, and measures the frequency of the Rabi oscillations, with 
quantum or classical character. This provides a unique opportunity to investigate a basic phenomenon of quantum 
electrodynamics in a condensed matter experiment. By changing the external magnetic field applied to graphene, 
the coupling between the atom and electromagnetic field in the JC model can be tuned continuously, facilitating the 
exploration of various regimes, the quantum to classical crossover. 

In quantum optics, one has the freedom to prepare the initial state of both the atom and the electromagnetic 
field^iiiS. The atom is usually prepared in its excited state, and the field is prepared in a number state or in a coherent 
statei^i^^. Then one can study the time evolution of the atomic population, which exhibits Rabi oscillations, when 
jumping between the ground and excited state, causing collapse and revival phenomenon. However, qualitatively 
different behaviour describes chaotic or thermal fields^. Quiescent periods and interfering revivals are also present. 
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FIG. 2: (Color online) The real time evolution of Cxx{t) is shown at T = 0, taking both valley and spin degeneracies into 
account. We introduced a cutoff D, and the number of levels is measured as D = V\/N, corresponding to different magnetic 
field strengths. The left/right panels show iV = 10000 (blue)/iV = 100 (red) with D = VVN for {ij,,A)/V\/N = (0,0) [a/e], 
(0.4,0) [b/f], (0.4,0.2) [c/g], and (0,0.2) [d/h]. These structures correspond to thermal field induced random oscillation in 
quantum optica^. 




FIG. 3: (Color online) The real time evolution of Cxx{t) is shown for long times for T = 0, (a*. A) — (0,0). We introduced 
a cutoff D, and the number of levels is measured as D = Vy/N, corresponding to different magnetic field strengths with 
= 10000 (blue) and A'^ = 100 (red). Collapse and revival shows up with time similarly to the thermal field Jaynes-Cummings 
model. The revivals gradually get wider and overlap. The presence of thermal revivals are related to the finite average boson 
number in the Jaynes-Cummings model, which translate to a finite cutoff in the Dirac case. As opposed to Fig. [5J these revivals 
at long times are caused by the finite cutoff. 

but the resulting pattern of oscillations follows an apparently random evolution. A coherent state, which is the 
most classical single mode quantum state, is strongly peaked around the average boson number n. Therefore, Rabi 
oscillations mainly involve frequencies around Efia, and the collapse time is independent of the field strength, i.e. h. 
As opposed to this, a thermal or chaotic field has a broad, monotonically decreasing distribution of boson states, and 
can be represented as a mixture of coherent states with a Gaussian distribution of the mean values. As a result, the 
very wide range of boson numbers gives such a broad distribution of Rabi frequencies that almost no trace of coherent 
oscillations remains after ensemble averaging, and the resulting initial collapse time depends on the field strength h. 
The revival times depend on n for both initial conditions. 

For Graphenc, an arbitrary preparation of the initial states is not accessible, but requires thermal, ensemble av- 
eraging. In this respect, it is closer to the second type of thermal initial condition for the JC model. The average 
boson number in the JC model corresponds to the total number of fermions in graphene, determined by the chemical 
potential and the cutoff. This can be introduced by the energy scale D = V\/N + 1, above which we neglect all states 
(with 71 > N). We mention that the inclusion of a cutoff is required to obtain correctly the f-sum rule for graphene»2i. 

Using this prescription, we investigate the real time evolution of the longitudinal current-current correlation function 
based on Eq. The results are shown in Fig. [21 including valley and spin degeneracies. Similarly to observations in 
quantum optics^ii^, the initial collapse is followed by a revival of oscillation, which are also sensitive to the presence 
of finite fi and A. They both enlarge the quiescent period after the short time collapse, and cause additional step-like 
structures in the envelope of oscillations with a period of max(/.t, A)Tr/vpeB. For longer times, collapse and revival is 
observable in Fig. [3l which gradually become wider and overlap. This revival time depends on the value of the cutoff 
like 2tt\/ N + 1/V = irD/vpeB, as is apparent from the figure, and is controllable by the magnetic field. 

The Hall response evolves differently. In the DC limit, by taking both valley and spin degeneracies into account, it 
produces the unconventional quantum Hall steps as a function of fi. In connection with the JC model, these steps. 
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FIG. 4: (Color online) The real time evolution of Cxy{t) is shown for T = 0. The explicit value of the number of levels (A'^) 
does not influence the resulting pattern. The chemical potential varied as fi/V = 1.2 [a/d], 3.2 [b/e] and 10 [c/f] with A = 
(left panel, blue) and A — 2V (right panel, red), and the frequency of the envelope function is v%eB/fi, the cyclotron frequency 
of massless Dirac fermions. For T = and A > ^, ImCxy{t) = 0. Note the different horizontal scales! 



occurring when ^ coincides with a Landau level energy, correspond to the bare Rabi frequencies, which can be revealed 
in a static, DC experiment, but still gaining information about the dynamics of the system. The microwave Hall 
conductivity obeys to the same selection rules and possesses the same transition matrix elements as the longitudinal 
conductivity. However, the extra phase factor in Eq. ([2|) is responsible for a different behaviour. 

The time evolution of the Hall correlator is shown in Fig. [H and turns out to be independent of the applied 
cutoff scheme, hence being universal. It exhibits coherent Rabi oscillations, which vanish at the Dirac point (this 
is equivalent to the statement, that the Hall conductivity is zero exactly at the Dirac point). Upon increasing /i, 
oscillations shows up with beating property, observed in quantum optics as well^. The characteristic frequency of the 
envelope of the oscillations for /i ^ V, A is V^/2fi = v^eB/fi, which is the cyclotron frequency of massless Dirac 
fermions^. Even though the spectrum is linear, the finite chemical potential provides us with an energy scale for the 
cyclotron mass. Only the energy levels with \Ena\ < M contribute to the Hall response at T = 0. Therefore, for 
A > fi, hnCxyit) = and no oscillations are present. This means, that only a finite, narrow range of frequencies 
determine the Rabi oscillations, similarly to the coherent field case in quantum optics, which lead to beating and 
non-chaotic collapse and revival. Thus, although the initial field state is always a thermal one in graphene, both 
coherent and chaotic Rabi oscillation can be observed in different quantities, together with collapse and revival. Thus, 
the correlation function C{ip, t) measures in principle the crossover from thermal to coherent behaviour with changing 
(f, and provides us with the unique opportunity to observe the crossover by not changing the initial field, but by 
measuring a different component of the current. 

Time-resolved current- voltage measurements can reveal the real-time dependence of correlation functions, similarly 
to the current-voltage characteristics of Josephson junctionsS^. By applying a sharp current or electric field pulse 
as E = EoSit), the resulting current through the sample parallel or perpendicular to Eq after the initial pulse 
is directly related to the above correlation functions. More precisely, within linear response theory, it follows as 
(jait)) = —2e^VpEQ Jq dt'lmCxa{t') with a = x or y. This opens up the possibility to reconstruct the time-dependence 
of ImCa;a(t'), or to obtain the frequency dependent longitudinal and Hall responses. Another method invokes the 
fcmto/atto-second laser pulse technique. After shining the sample with a short laser pulse, the measurement of 
transmittance or reflectance in time is determined by the C(t) correlation functions. The most conventional way to 
deduce these correlation functions is provided through optical conductivity or current fluctuation measurement. Via 
the fluctuation dissipation theorem, they contain the same information, and upon Fourier transforming from frequency 
space to get the real time dependence, one is expected to be able to observe the presence of thermal Rabi oscillations. 
Similar measurements have already been carried out without magnetic fiel d^^i^^ , by exploiting the tunability of the 
carrier concentration with gate voltage. 

The above calculations can easily be extended to other correlation functions such as thermal conductivities. We also 
may speculate that various extensions to single layer graphene such as bi- and multilayer structures^ (with spin-orbit 
coupling) can be mapped onto different multimode, multiatom and non-linear versions of the JC model. 

We have shown that the equivalence of the Hamiltonians of graphene in magnetic field and of the JC model influences 
their correlation functions as well, causing both thermal and coherent Rabi oscillation in the electric response of 
graphene. Finally we speculate that Rabi oscillations and Zitterbewegung are two closely related phenomena named 
differently in different fields of physics, both arising from the coupling of positive and negative energy states. 
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